test analytical predictions [17]. 
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Figure 5: Coefficient of the log 2 (6 2 ) versus the temperature for different values of the A 
parameter from a two parameter fit. The same symbols as in figure 3. 

4 Conclusions 

This comprehensive study reveals clear signals of a transition both in the static and dy- 
namic properties of the system. The transition seems to occur at the same temperature 
whether determined by static or dynamic methods. The transition temperature does how- 
ever depend on A, being close to k/tt only at small A. We find the general behavior 
predicted by RG arguments, that is: linear behavior of the dynamic exponent, z oc r, and 
quadratic behavior of the coefficient of log 2 , 6 2 oc r 2 . Detailed numerical agreement of the 
coefficients, especially in the case of the quadratic behavior, is lacking (by a factor of ~ f 4) 
while in the dynamic case the agreement is reasonable. Finally, an explanation of the A 
dependence of all the quantities measured is needed. 

The essential difference of opinion on the low temperature phase is between analyses 
based on the RG that predict a log 2 form for the correlator, and analysis based on a varia- 
tional approach that yield a log. The RG approach has been criticised because the solutions 
are unstable with respect to breaking of replica symmetry [12], however the significance 
of this is not clear since the techniques for dealing with replicas are not sufficiently devel- 
oped to deal effectively with non mean held situations. On the other hand, the variational 
approach has only been calculated for a gaussian ansatz which corresponds to leading or- 
der in some 1/N expansion. Recent results [16] using RG arguments for an N-component 
version of the Random Phase Sine Gordon model are interesting in that they calculate the 
coefficient of the log 2 term to be order 1/N 3 . We hope that further analytic work in both 
approaches to understand the smallness of the log 2 coefficient seen in this and all other 
simulations, will lead to a resolution of the puzzle. 

Finally we believe that we have shown that it is possible obtain reliable and accurate 
numerical data for this rather contentious subject that can be extended and used to further 
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Figure 4: Coefficient of the log (&i) versus the temperature for different values of the A 
parameter from a two parameter fit. The same symbols as in figure 3. The error bars are 
smaller that the size of the points. 

jack-knife. Figures 4 and 5 show the resulting values of the coefficients for A = 0.5 and 2.0, 
and also the data for A = oo from reference [8]. 

The hrst coefficient bi continues to follow its high temperature behavior down to low 
temperatures showing only small deviations from linearity. The second coefficient b 2 is 
sensitive to the transition and grows at low temperature, it does not however display the 
same kind of linearity evident in the dynamic case. It is therefore slightly more difficult to 
determine the critical temperature though it is clear that this is A dependent. 

The form (8) is the result of RG calculations [I] which predict a universal behavior 
for the second coefficient b 2 = 8r 2 + 0(t 3 ) [2]. This quadratic dependence on reduced 
temperature can be tested, and we show in Table I (right) the results of a fit 



where we have only used points for which b 2 is further than twice the error away from zero. 
The fit is good, but in this case the value of the coefficient is not even close to the RG 
prediction. 

The transition identified in the dynamic and static measurements should occur at the 
same temperature. This follows from the RG analysis and also from a variational calcu- 
lation of the entropy of metastable states [14, 15]. We see from Table I that both static 
and dynamic measurements give similar critical temperatures, and most importantly, that 
these shift simultaneusly with different values of A. 




(9) 
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Dynamics 


Statics 


A 


Interval 


a D 


T D 


Interval 


a s 


T s 


0.5 


[0.4,0.5] 


2.5(5) 


0.53(2) 


[0.4,0.5] 


0.69(12) 


0.58(2) 


[0.35,0.5] 


3.1(3) 


0.52(1) 


[0.35,0.5] 


0.67(7) 


0.58(1) 


2.0 


[0.5,0.75] 


3.8(2) 


0.78(1) 


[0.5,0.75] 


0.54(6) 


0.77(2) 


[0.6, 0.7] 


4.3(6) 


0.76(2) 


[0.55,0.75] 


0.44(12) 


0.78(4) 


oo 


[0.55,0.8] 


2.7(3) 


0.85(3) 


[0.5,0.80] 


0.44(7) 


0.87(2) 


[0.6, 0.75] 


2.8(6) 


0.83(4) 


[0.6,0.80] 


0.35(5) 


0.90(2) 



Table 1: Results of linear fit to z(T) (dynamics) and quadratic fit to b 2 (T) (statics) for 
three values of A and for different sets of temperatures (intervals). See text for notation. 
All the fits in this table have x 2 /d.o.f 1. 

We have shown in this section that in the vicinity of the critical temperature the z- 
exponent has a linear dependence on reduced temperature but that the critical temperature 
depends on A. The slope is in tolerable agreement with the predicition from RG. 



3 Statics 

In this section we discuss the equilibrium form of the correlation function along the lines 
of previous work on the discrete gaussian model [8] and shall repeat some of the arguments 
used there. The main priority to to ensure that thermal equilibrium is actually obtained. 
An annealing scheme is used and it is important to include the large Monte-Carlo moves 
discussed above, even so, very long times are needed. Typically for each temperature, 
we thermalise for 10 5 sweeps, measure over a period of 2.10 4 sweeps, and repeat this 
cycle several times. The system is taken to be thermalised when we reach a situation in 
which subsequent measuring cycles show no systematic drift and fluctuations within each 
measurement cycle are similar to those between cycles. 

At high temperatures the correlator is accurately fitted by a single log, or more precisely 
by the lattice version: 

i L-1 L-1 i _ /27rni| \ -, 

U 2^^ 1 ^ 2-cos(^)-cos(^)-2 7 r gl 2^ J ' 10 

where the symbol ~ holds for L ^> 1. The coefficient follows the temperature to within 
1%, in agreement with all theoretical predictions [4, 11]. 

As the temperature is reduced below some value this fit becomes substantially worse. 
The degradation cannot be ascribed to short distance lattice effects since a similar worsen- 
ing occurs for fits in which short distance points are omitted. We emphasise that the long 
distance finite size effects are well under control in equation (7). Our results are not sen- 
sitive enough to determine the functional form of the correction necessary but theoretical 
prejudice suggests a fit of the form [1]: 

C{r) = b 1 P L {r) + b 2 P L {rf. (8) 

This fit works very well for all lambda, improving as equilibrium is approached, and with 
a value of % 2 that is small and does not vary significantly throughout the temperature 
range. The fit is made by an exact minimisation procedure and errors are determined by 
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Figure 3: The dynamical critical exponent, z, against temperature for A = 0.5, 2.0 and oo. 
The gaussian value is marked with a horizontal line. We use squares and a dotted line for 
A = 0.5, diamonds and a dashed line for A = 2.0 and crosses and a continuous line for 
A = oo. The lines are only to guide the eye. 
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Figure 2: Growth of correlation length with time (A = 2.0, T = 0.5). 



explicitly by comparing with the slower dynamics based on a simple move taken from a 
single flat distribution centred around zero, and in fact have used this alternative move 
for our determinations of z. The dependence of z on temperature is shown in figure 3 
for A = 0.5,2.0 and oo. At high temperatures the dynamic exponent takes its gaussian 
value of 2, marked by a line in the figure. Below, but in the vicinity of the transition the 
deviation from the gaussian value appears to be linear in the reduced temperature. This 
observation would be in accord with the RG prediction for T < Tc [10] 



where the formula is valid for small A with Tc = 2/tt and 7 is Euler's constant. If the 
linearity continues to hold for large values of A we can determine transition temperatures 
reasonably accurately and it is clear that these depend on A. In any event, the large A 
data is certainly not compatible with Tc = 2/tt. The formula (5) yields a A independent 
slope with numerical value 3.56. We can try to extract the slope and critical temperature 
assuming a fit of the form 



The results for ap and Tp are reported in Table 1 (left). Different sets of points have been 
used to demonstrate that the fit is stable, but we always ignore points within twice their 
error of the gaussian value (z = 2). 

An alternative type of comparison is available with results for the dynamical exponent 
based on the response to a small driving force F. In the low temperature phase RG 
arguments give the relation between applied force and resulting velocity as [3]: (dt<f>) ~ 
t 2z F 1 + 2zT 5 where r is the reduced temperature. It is reassuring that our data are in good 
numerical agreement, both at small [5] and large [9] A, with such a different approach based 
on Langevin dynamics. 




(5) 




(6) 
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Figure 1: Correlation functions, C t (r) } shown at equal time intervals (bottom to top: 
t =2000, 4000, 6000 and 8000 sweeps), and the equilibrium correlator, C asy (r), as the 
continuous line (A = 2.0, T = 0.5). 

where ( • ) and ~ denote thermal and disorder averages respectively. 

2 Dynamics 

It is convenient to plot the correlator against log r as in figure I, the top curve represents 
the equilibrium form. This equilibrium curve is fairly linear except for the finite size effects 
and a small component of a different large distance behavior which will be the subject of the 
next section. The other curves in figure I illustrate the approach to this final equilibrium 
form, the short distance part of the correlators follows the equilibrium curve but a plateau 
is reached beyond a distance that can intuitively be identified as the correlation length £(t) 
at that time. 

This intuitive definition of the time dependent correlation length is numerically robust 
since the plateaus are very well defined. Knowing the asymptotic form of the correlator 
(C asy (r)), it is simple to obtain the correlation length £(t) from C asy (£(£)) = Plateau[Cf(r)], 
where C t (r) is the correlation function at the time t. Starting from a constant configuration 
(i.e. (pi = Const.), we measure how the correlation length grows with time. Provided £(t) 
does not approach the system size, we find that dynamical scaling is well obeyed: 

m ~ f xiz . (4) 

We show this growth in figure 2 on a log-log scale. The errors on each point arise 
from variations between different samples. We emphasise that this scaling behavior is 
independent of the detailed choice of local Monte-Carlo dynamics. We have checked this 
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1 Introduction 



The Random Phase Sine-Gordon model has various physical interpretations of interest 
such as a crystalline surface with disordered substrate [1] and an array of flux lines in 
a superconducting film subject to random pinning and parallel magnetic field [2] . The 
model has recently been the scene of considerable activity, both theoretical and numerical 
(see for instance [3, 4, 5, 6, 7, 8, 9]), and in neither area is there agreement on the details 
of the low temperature phase. In this letter we present a comprehensive numerical study 
of both dynamic and static aspects of the model for a range of values of the interaction 
parameter A. We also compare with the disordered discrete gaussian model which was the 
subject of earlier work [8]. This study reveals simultaneous movement of the static and 
dynamic critical temperatures as A is varied, a phenomenon not previously visible in more 
restricted studies. 

The Hamiltonian of the Random Phase Sine-Gordon model is: 

# = ^) 2 - A5>s(27r(& - rn)), (1) 

<ij> i 

where < ij > denotes nearest-neighbors, cf>i is a continuous variable and rji £ [0, 1) is the 
quenched random disorder. 

The discrete gaussian model with a disordered substrate is related to the A — > oo limit 
of this model. 

# = h)\ (2) 

where the field is now integer valued, up to shifts by the disorder: cf>i = rii + rji ;rii £ Z. 

We employ a Metropolis Monte-Carlo algorithm, and for the discrete gaussian model 
the proposed change in the field is naturally ±1. For the continuous model the move is 
chosen from a distribution consisting of three flat regions, each of width e and respectively 
centred about and ±1. This form of move is far more efficient than one taken from a 
single flat distribution centred at zero, and makes the connection between the two models 
clearer. We have fixed the relative height of each region so that about half the moves are 
large. The widths e are adjusted at each temperature and value of A to give an acceptance 
rate of approximately 50%. 

We have used the computer APE [13], simulating in parallel 256 systems composed 
of 128 pairs with distinct disorder (samples). Each pair of uncoupled replicas with the 
same disorder was subject to different thermal noise, and in the standard procedure for 
disordered systems the overlap between the copies was monitored. The disorder average 
was taken over all 256 systems. 

The size of the system is always 64 2 , studies of the discrete gaussian model for larger 
sizes were reported in [8]. We fix n = 2 so the critical temperature according to theory 
valid for small A is Tc = 2/tt. Various values of A have been considered up to A = 3.0, 
thermalisation becoming increasingly difficult as A increases. Here we present results for 
A = 0.5,2.0. 

In all cases the energy relaxes almost immediately to its final value and we therefore 
conclude that local equilibrium is achieved quickly, leaving the long time dynamics to be 
that of domain reajustment. A direct visual inspection of the configurations confirms this 
but a detailed investigation of the process requires the correlation function defined by: 

C(r) = < (<f> r - <^o) 2 >, (3) 
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Abstract 

We have performed comprehensive numerical simulations of the Random Phase 
Sine Gordon Model, studying both statics and dynamics for various values of the 
coupling. The glass transition can be seen both in static and dynamic signals at a 
temperature that depends on the coupling. Our results agree qualitatively (statics) 
and quantitatively (dynamics) with Renormalisation Group predictions. 
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